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INVESTIGATION OF THE STARTING CHARACTERISTICS OF THE 
M-l ROCKET ENGINE USING THE ANALOG COMPUTER 
by John R. Szuch, Leon M. Wenzel, and Robert J. Baumbick 
Lewis Research Center 

SUMMARY 

The purpose of this investigation was to study the starting characteristics of the M- 1 
engine system by using an analog computer simulation. It is further intended that this 
report serve as a guide to others undertaking liquid rocket propulsion system studies on 
the analog computer. The formulation of the M-l system equations is discussed in de- 
tail, and a summary of the resultant equations is included in an appendix. To evaluate 
the starting characteristics of the engine, the times required to reach 90 percent of the 
rated thrust chamber pressure and fuel pump stall margin were selected as figures of 
merit. The following information was obtained from the analog simulation: 

(1) The effect of ambient pressure on the engine start transient was negligible with 
pressures ranging from zero to atmospheric. 

(2) Stall margin could, in general, be increased by increasing the recirculation valve 
area at the cost of increased acceleration time. Stall margin was found to be a slight in- 
verse function of recirculation line diameter. A decrease in acceleration time, together 
with an increase in stall margin, was obtained by using a 4-inch line with a recirculation 
valve area and spring rate greater than nominal. 

(3) Failure of the recirculation valve to close during the start transient resulted in a 
steady- state thrust chamber mixture ratio near stoichiometric. 

(4) For any given inlet line configuration, thrust chamber to propellant tank pressure 
ratios can be obtained from engine impedance characteristics measured on the analog. 
Over the expected range of possible POGO-stick-type oscillations (0 to 20 cps), the en- 
gine impedances were approximately constant with the dominant dynamics expected from 
the inlet lines. 


INTRODUCTION 


The purpose of this investigation was to study and predict the dynamic characteris- 



tics of the M-l engine system by using an analog computer simulation. It is further in- 
tended that this report serve as a guide to others undertaking studies of liquid rocket pro- 
pulsion systems on analog computers. 

A computer simulation, when properly used, can be a powerful tool in guiding an en- 
gine development program. It provides an easy and economical means for evaluating 
various design approaches and forewarns the designer of possible problem areas before 
costly hardware is developed and subsequently scrapped. Once a qualitative design is es- 
tablished, the system may be "tuned” for optimum performance by varying the system 
parameters about their design values. 

Although the digital computer provides the design engineer with higher numerical ac- 
curacy, the analog computer offers the following advantages over the digital as a design 
tool: (1) the system equations are solved continuously rather than at discrete intervals of 
time; (2) the engineer is able to observe the effects of changes in the system as they oc- 
cur rather than in numerical form at some later date; and (3) visual displays of system 
behavior, such as engine start transients, etc. , give the engineer an added "feel” for 
his problem. 

Since a major portion of the M-l engine hardware was still in the design stage, de- 
sign information, rather than experimental data, was used to simulate the system compo- 
nents. Wherever possible, information from other engines was used to aid in the simula- 
tion. For example, the simulation of the RL-10 engine, for which there was excellent 
agreement with test data (ref. 1), provided invaluable simulation experience and informa- 
tion. Several of the approximations and assumptions, described in the analysis, were 
based on experience gained through the RL-10 engine simulation. Acknowledgement is 
made to Aerojet General Corporation for their cooperation in providing the necessary de- 
sign information in addition to results from their digital computer simulation of the M-l 
system. 

Under limitations of using design data rather than test data, the work reported here- 
in pertains to engine simulation in general and the development of the M-l simulation in 
particular. To aid in the design of facilities for engine testing, the effect of ambient 
pressure on the start transient was determined for pressures ranging from zero to atmo- 
spheric. To optimize the design of the fuel pump recirculation system, parametric 
studies on line diameter, valve diameter, and valve spring rate were conducted. In ad- 
dition, a failure analysis was performed on the recirculation valve to determine the criti- 
cality of such a failure. Engine transfer functions were determined that will be useful in 
vehicle design studies once tank and inlet line data become available. 


EQUIPMENT AND PROCEDURE 

A schematic of the M-l system is shown in figure 1. The engine is powered by liquid 
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Figure L - Schematic of M-l rocket engine system. 


hydrogen and liquid oxygen supplied by the two pumps and is regeneratively cooled by the 
liquid hydrogen that passes through the cooling jacket before entering the thrust chamber. 
Initial pump spinup is accomplished by supplying high pressure helium gas from the start 
bottle to the turbines. After the gas generator ignites, the helium flow rate gradually de- 
creases as the gas generator burnt products supply the energy for acceleration and steady- 
state operation of the turbines. Turbine discharge gases are expelled through the nozzle 
extension skirt to provide skirt cooling. Rising fuel pump discharge pressure closes the 
fuel pump recirculation valve, which is initially held open by spring forces. This valve 
provides both prestart pump cooling and a low flow impedance downstream of the fuel 
pump during the initial phase of the start transient. 

The M-l analog simulation utilized the analog computer facility at the Lewis Re- 
search Center. Three computer consoles were used with all mode controls slaved to one 
console. Approximately 250 amplifiers, 50 multipliers, 20 variable diode- function gen- 
erators, and 5 fixed diode- function generators were used in the M-l simulation. Start 
transients were recorded on three 8-channel strip-chart recorders. Transient traces 
were plotted on both oxidizer and fuel pump maps when needed by using X-Y plotters. 
Program expansion and refinement, whenever necessary, was accomplished by using an 
additional console. 

Since the computers used are 100-volt devices, all system variables were scaled so 
as not to exceed that value. Included in appendix B is a summary of the voltage scaling 
used in this simulation. To facilitate the programing and recording of the start transient, 
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the system was slowed down (time scaled) by a factor of 100. The sequencing of system 
valves was accomplished by using pump discharge pressure to drive comparator-driven 
relays, while the start valve opening and closing was accomplished by manually throwing 
a switch on the computer control panel. 

Because of operational problems, the start transient could not be initiated on the 
computer from an absolute zero condition. For example, an attempt to compute initial 
pressure ratios would result in fallacious results since the electronic dividers are inac- 
curate for very small input voltages, small errors being greatly magnified. To avoid 
this situation, the problem was started with a small helium flow (0. 5 percent of maximum) 
through the hot gas ducts with appropriate initial conditions on the problem variables. 

The engine transfer functions were obtained experimentally by using the computer 
model, minus the inlet lines, and a sinusoidal transfer function analyzer. The latter de- 
vice generates its own test signal (input) and receives whatever output signal is of interest 
(e. g. , thrust chamber pressure). The analyzer then computes the magnitudes and phases 
of the Fourier transforms of both input and output signals at test frequencies selected by 
the operator. Since the problem was time scaled by a factor of 100, and since the fre- 
quency range of interest was 1 to 100 cps, the analyzer was operated from 0. 01 to 
1. 0 cps. 


ANALYSIS 

The first step in developing any analog program is the formulation of an analytical 
model. This formulation will, in general, result in a set of simultaneous equations. 

Since the analog computer can handle only ordinary differential equations, certain ap- 
proximations and manipulations are necessary when preparing a problem for the analog 
computer. For example, one -dimensional heat flow problems result in partial differen- 
tial equations, which can be handled by means of finite- difference techniques (ref. 2). 

In a complex system simulation, such as the M-l, approximations must also be made to 
conserve computing equipment. Engineering judgment is required to economize on equip- 
ment where the result will not invalidate or seriously compromise the overall problem 
accuracy. 

A block diagram of the M-l system is shown in figure 2. To facilitate the program 
development, the M-l engine system was considered as being a composite of three sec- 
tions. These sections were the helium and hot gas driving system, the main propellant 
feed system, and the cooling jacket and main engine assembly. After deriving the neces- 
sary equations, each section was built up on the computer independently. The simulation 
was then completed by closing the various interconnecting loops. In the following para- 
graphs, these sections will be discussed in detail in the order just mentioned. 
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Figure 2. - 



diagram of M-l engine system. 





































Helium and Hot Gas Drive Section 


On startup, the M-l turbines derive their power successively from the flow of helium 
from the start bottle, a mixture of helium and gas generator burnt products, and finally, 
all burnt products. The burnt products are supplied by the gas generator, which is fed by 
liquid hydrogen and oxygen tapped from the main propellant feed lines. All losses be- 
tween the main-line tapoff and the gas generator chamber were lumped together because 
of equipment limitations. These losses (lines, valving, manifolds) were assumed to be 
proportional to the square of flow rate (ref. 3). 

The gas generator valve configuration is shown in figure 3. Both the fuel and oxidant 
valves are driven by the same actuating piston and are opened by the rising fuel pump dis- 
charge pressure. An economical simulation of the valve dynamics was made by neglect- 
ing valve forces, assuming that the timing orifice is always choked, and neglecting the 


change in the piston actuating volume V PA 


due to the piston travel. 


These assumptions, 


together with the prescribed relation between valve contours and valve stroke, led to the 


following equations, which are easily implemented on the computer: 


Valve stroke: 

y GGV = K 1 P dFP dt ' *2 ^ 

Gas generator oxidant flow rate: 

/ \ r / x i 1 / 2 

W GGO = ( C d A ) 0 L 2 ^o( AP GGv) 0 J < 2 > 

Gas generator fuel flow rate: 

. r -.1/2 

W GGF = ( C d A ) F [_ 2gp F( AP GGv) F J 

All symbols are defined in appendix A. 

The continuity equation (ref. 3) was used to compute the gas generator chamber pres- 
sure from the net instantaneous weight flow into the chamber, assuming a perfect gas 
after ignition and no combustion dead time. The gas properties in the gas generator 
chamber (i. e. , C*, R, and y) were computed from the percent fuel entering the gas gen- 
erator by using curve fitting techniques. Gas generator ignition was delayed until the 
liquid oxygen manifold was filled. The filling time was determined by integrating the gas 
generator oxidant flow rate until the weight stored in the manifold volume was equal to 
the capacity of that volume (Po^man^ • 

An engine balance orifice (fig. 2) is located between the gas generator and the fuel 
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turbine manifold. The sizing of this orifice affects the distribution of pressures across 
the fuel and oxidant turbines, which allows the thrust level to be adjusted by means of the 
pump driving powers. The orifice was assumed never to choke, and the flow rate through 
the orifice was computed by using equation (4). To account for the backflow of helium 
through the orifice during the first portion of the start transient, the fluid density at the 
orifice was computed from the average of the upstream and downstream pressures and 
the temperature downstream of the orifice (i. e. , at the helium - O/F mixing region): 

^ORF = ( C D A ) ORF ( 2gp ORF AP ORf) ^ 

For the computation of both turbine flow rates, equation (4) proved to be inadequate due 
to the higher flow velocities. In the case of a high Mach number (M > 0. 3) restriction 
(ref. 4, pp. 83 to 89), the flow must be treated as compressible, which requires a know- 
ledge of static downstream pressures. For both turbines, it was necessary to relate 
static pressure to total pressure at the turbine exit. It was assumed that the static to 
total pressure ratios at the turbine exits were constant: 


= constant = K r (5) 

r 

over the range of Mach numbers encountered during the start transient (M < 1.0). This 
assumption allowed the computation of the turbine pressure ratios as follows: 


fd _fd 

P d' P d 



( 6 ) 


The use of the energy equation led to the following expression for the flow rate 
through the turbines in terms of the pressure ratios (ref. 4) at the nominal mixture ratio 
in the gas generator: 


W, 


TUR 


P u A TUR 

<% T / /2 




( 7 ) 


Equation (7) implies a pure impulse machine with no static pressure loss across the 
rotating elements. Although equipment limitations prevented any modification of the tur- 
bine flow rates due to variations from this condition, the resulting turbine torques were 
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continuously adjusted through turbine efficiencies, which were computed from equations 
supplied by the contractor. 

To account for variations in the turbine flow rates due to changes in fluid properties, 
a curve fitting technique was used to obtain a first-order correction factor: 

F <V M > - k 3 + k 4 < 8 > 

y M 


Combining equations (7) and (8) resulted in the following expression for the flow rate 
through the turbine (either fuel or oxidizer turbine) : 


W 


TUR 


F(y M )^TUR 


r 


(9) 


The flow rate through the propellant utilization valve (see fig. 2, p. 5) was computed in a 
similar fashion. This valve controls the mixture ratio in the thrust chamber by control- 
ling the flow rate through the fuel turbine. The temperature drop across the propellant 
utilization valve was considered negligible. 

In order to compute the main line flows, pump speeds, etc. , it was necessary to 
compute the turbine torques. The following equations, derived from energy considera- 
tions and curve fitting techniques, summarize the methods used in computing the turbine 
torques and downstream temperatures (ref. 4): 


Ah ID " c pM T u 


Ah 


1- |K — - 

P u 


ACT “ ^TUR Ah ID 


M 


Itur - K s(2) * K e(“) 


u = M rTURN 


c » (2gJ ih ACT ) 1/2 
Ah 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


ACT 


L - W 

xUR 2v W TUR 


N 


(15) 



c pM T d “ c pM T u " Ah ACT 


(16) 


Fuel turbine: 


Oxidizer turbine: 


c pM T u ~ c pM T M 


c pM T u “ c pM T M 


^TUR Ah ACT \ 
^TUR + ^PUV/ F 




( L TUR " L P^ dt 


(17a) 


(17b) 


(18) 


The determination of turbine torques, therefore, required continuous knowledge of 
the helium - O/F mixtures* temperature and. gas constants. These were computed by 
averaging (by weight) the contributions of the helium and burnt products. The following 
equations summarize this averaging technique which is developed in detail in appendix C. 


0 _ W GGO + W GGF 

(19) 

^GGO + ^GGF + ^HE 


C pM = ^ c pBP + ^ " ^ C pHE 

(20) 

h M = ^ h BP + ^ “ ^ h HE 

(21) 

r m = ^^p + (- 1 “ ^ r he 

(22) 


The nozzle extension skirt was considered to be a choked orifice with the flow rate 
through the skirt computed from the following equation: 


W 


A SK P dOT 


SK 


( R M T dOT^ 


1/2 


- 2 \^M-^/ 2gyM \^ 

+ V VM + V 


(23) 
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The helium gas is supplied by a start bottle (see fig. 2) initially at pressure and 
temperature T^. Upon opening of the start valve, the helium gas was assumed to expand 
isentropically. The pressure and temperature were then related (ref. 5) through the fol- 
lowing equation: 


T p (1 ' y HE )/r HE _ T p (1_r HE )/y HE (24 , 

T HE*HE ~ T i P i W 

The start valve was considered to be variable area restriction with the flow rate com- 
puted from equation (7). The static downstream pressure was assumed to be equal to the 
total pressure downstream of the valve (M « 1. 0). The dynamics of the line between the 
check valve and the fuel turbine manifold were neglected. All losses associated with the 
start valve and check valve were lumped together and assumed proportional to (W RE ) . 


Main Propellant Feed Section 

Finite difference techniques were used to simulate the suction and discharge lines to 
the main engine with representative line lengths chosen for the suction lines. The suction 
lines (both fuel and liquid oxygen) were represented by symmetrical T-shaped sections 
and the discharge lines by rr-shaped sections. The T-shaped representation of the suction 
lines allowed explicit expressions to be written for the pump flow rates. 

To avoid fuel pump stall during startup and to provide prestart pump cooling, a re- 
circulation valve is located between the pump discharge and fuel tank. Because of the de- 
pendence of stall margin and acceleration time on this valve and its associated line, a 
careful simulation was made. The line was represented by four 7r-shaped sections, util- 
izing the aforementioned finite difference techniques. The flow rate through the valve 
was computed from the following equation: 

^RCV = ( C d A ) RCV ( 2gp F AP RCv) 1//2 ( 2 5 ) 

The fuel pump recirculation valve is spring loaded open and is closed by increased pump 
discharge pressure as the engine bootstraps. 

The thrust chamber valves were similarly treated, when it was assumed that liquid 
conditions existed at the valves. These valves are initially held closed by spring forces. 
The opening rate of each valve is a function of its respective pump discharge pressure. 
Both thrust chamber valves are rigidly connected so that in effect the opening rate of both 
valves is a function of the average of the two pump discharge pressures (see eqs. (D51) 
to (D56) in appendix D) . 
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Figure 4. - Estimated pressure rise char- 
acteristics for eight-stage axial fuel pump. 


Figure 5. - Estimated torque characteristics 
for eight-stage axial fuel pump. 


The pump characteristics (i. e. , pressure rise against pump flow rate and torque 
against flow rate) supplied by the contractor were represented by curve fits (eqs. (26) 
to (31)). Pump efficiencies were included in the expressions for pressure rise and torque. 
The oxidizer pump is a single-stage centrifugal type and the fuel pump is an eight- stage 
axial design. The fuel pump was represented by two sets of equations. One set was valid 
for the unstalled region while the other represented the pump while in stall. Only one set 
of equations was used to represent the oxidizer pump since a large stall margin was ex- 
pected with the present design. The fuel pump pressure rise and torque characteristics 
are shown in figures 4 and 5. Transitions between the stalled and unstalled regions were 
made possible by continuously comparing W/N for the fuel pump to the critical value of 
W/N, which is characterized by the stall line (figs. 4 and 5). Because of lack of experi- 
mental data, the stall line is not defined for speeds less than about 2 5 percent of nominal. 
The following equations were used to compute head rise and torque for the two M-l pumps: 


Fuel pump: 
Unstalled: 


N 2 ^ 8n 



AP 

N 2 


= K 


10 


11 N 




(26) 

(27) 
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Stalled: 



( 28 ) 

(29) 

(30) 

(31) 


The oxidizer pump pressure rise and torque characteristics are presented in figures 6 
and 7. 

The pressure drop from the thrust chamber valve on the fuel side to the jacket was 
assumed to be proportional to the square of the flow rate through the valve. The liquid 



0 50 100 150 200 

Percent- rated flow rate 


Figure 6. - Estimated pressure rise character- 
istics for single-stage centrifugal liquid oxygen 
pump. 



0 50 100 150 200 

Percent- rated flow rate 


Figure 7. - Estimated torque characteristics for 
single-stage centrifugal liquid oxygen pump. 
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oxygen side pressure drop was also assumed to be proportional to the square of the flow 
rate with the injector drop lumped with the line losses. 


Cooling Jacket and Main Engine Assembly 


The flow rate through the jacket was assumed to be compressible with heat transfer 
taking place between the combustion zone and coolant fluid. Based on RL-10 experience, 
the flow rate was assumed to be proportional to the pressure drop across the jacket 
(ref. 1). 

A one lump finite difference approximation to the jacket was used with the average 
temperature in the jacket computed from the net heat flow into the fluid. Under the pre- 
viously mentioned limitation on computing equipment, the following assumptions were 
made to facilitate the simulation of the cooling jacket (ref. 5) : (1) heat flow through the 
walls was proportional to the temperature difference across the wall; (2) heat flow through 
the coolant films was proportional to the eight-tenths power of the coolant flow rate; (3) 
the fuel injector flow rate was computed from equation (25); (4) the average temperature 
in the jacket Tj was equal to the average of the jacket inlet temperature and the fuel in- 
jector temperature; and (5) heat transfer between the environment and the outer wall was 
neglected. Figure 8 shows the one-dimensional model used to represent the cooling 
jacket. Heat transfer was considered between the combustion zone and coolant tube wall 



/ 

Coolant- / Combustion 
tube inner wall - zone 


and between the coolant and the coolant 
tube wall. Axial and circumferential heat 
transfer within the coolant tube walls, 
between the inner and outer sections, was 
neglected. Equations (D73) to (D89) were 
used to describe the heat-transfer process 
in the M- 1 cooling jacket. 



The thrust chamber pressure was 
computed from 


dP r 

T — + 

dt 




(32) 


Thrust chamber ignition was delayed until 
the liquid oxygen injector manifold and its 
feed line were filled. The filling time was 


Figure 8. - One-dimensional heat-transfer model for thrust chamber determined by using the same type of pro- 

coolant jacket 
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TABLE I. - RESULTS OF AMBIENT PRESSURE STUDY 


Ambient 

back 

pressure, 

P dsk’ 

psia 

Acceleration time 
to 90- percent-rated 
thrust chamber 
pressure, 

*90’ 

sec 

Thrust 

chamber 

ignition 

time, 

trc’ 

sec 

Gas generator 
ignition time, 

‘GG’ 

sec 

0 

1.49 

0. 69 

0.76 

.25 

1. 52 

.68 

.77 

1. 50 

1.49 

.69 

.76 

5.00 

1. 51 

.68 

.79 

14.70 

1.60 

.75 

.82 


cedure discussed in connection with the gas generator. 

Since this simulation was not intended as a study of combustion instability, combus- 
tion delays were not considered in the analysis of the thrust chamber. A summary of all 
the equations used in this simulation is given in appendix D. 


RESULTS AND DISCUSSION 

One of the fundamental objectives of the analog program was the determination of the 
sensitivity of the M-l start transient to variations in system parameters. To determine 
the requirements for an altitude start facility, the ambient pressure (back pressure on 
the extension skirt nozzle) was varied from zero to atmospheric. A trend toward longer 
start times was noted (table I). Since the data fell within a band of 0. 1 second, however, 
it was decided that it was not necessary to design an altitude start facility solely for the 
purpose of duplicating a start transient. 

The sensitivity of the engine start transient to the fuel pump recirculation configura- 
tion was also investigated on the analog computer. This study was made by varying para- 
metrically the recirculation valve area, the spring rate, and the recirculation line diam- 
eter and by measuring as performance criteria the time to 90 percent thrust chamber 
pressure and stall margin. This stall margin is illustrated in figure 9 in which a typical 
start transient is plotted on a fuel pump map. It is defined as the minimum value of 

. The selection of — as the figure of merit is consistent with nondimen- 
N 

sional procedures for evaluating pump performance, namely the measurement of Q/ANr, 
where Q is the volumetric flow rate through the pump. This method avoids the complex- 
ities involved in determining pump load lines over a transient range of operation. 

Figure 10 shows the effect of the recirculation valve area and line diameter on fuel 
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Time to 90 percent thrust chamber pressure, sec 



Figure 9. - Typical fuel pump behavior 
during M-l start transient 



Fractional percent of nominal recirculation valve area 

Figure 11. - Effect of recirculation system parameters on 
acceleration time. All recirculation valve spring rates. 



Fractional percent of nominal recirculation valve area 

Figure 10. - Effect of recirculation system parameters on 
fuel pump stall margin. Fixed recirculation spring rate. 

pump stall margin for a fixed (nominal) recir- 
culation valve spring rate. Although increas- 
ing the valve flow area resulted in a significant 
increase in stall margin (stall free accelera- 
tions were not possible with valve areas less 
than 0. 6 times nominal), the line diameter had 
an inverse effect that was possibly caused by 
the increased inductance of the smaller line. 
Figure 11 shows the sensitivity of the acceler- 
ation time to valve area for 4-, 6-, and 8-inch 
lines. Since the acceleration time is a mini- 
mum and relatively insensitive to valve areas 
for the 4-inch lines, this is obviously the best 
choice for rapid, stall-free operation. 

Figure 12(a) presents stall margin against 
recirculation valve flow area for a 4-inch line 
and various recirculation valve spring rates. 
Figures 12(b) and (c) present similar data for 
the 6- and 8-inch lines, respectively, fil- 
er easing the spring rate resulted in a signifi- 
cant increase in stall margin for a given line 
diameter, due to the sustained low flow im- 
pedance at the fuel pump discharge. However, 
the stiffer spring resulted in a decrease in the 
sensitivity of stall margin to line diameter as 
shown in figure 13. At 1. 2 times the nominal 
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Wpp /w F p\ 

Stall margin, -jj — l-n— I , (lb/sec)/rpm 



(a) Fixed recirculation line diam- (b) Fixed recirculation line diameter, 6 inches. (c) Fixed recirculation line diam- 
eter, 4 inches. eter, 8 inches. 


Figure 12. - Effect of recirculation system parameters on fuel pump stall margin. 



valve area 

Figure 13. * Effect of recirculation valve spring 
rate on stall margin sensitivity to line diameter. 


spring rate, stall margin is relatively insensi- 
tive to line diameter for all valve areas. Accel- 
eration time was found to be insensitive to the 
recirculation valve spring rate. Therefore, the 
acceleration time against valve area plot in fig- 
ure 11 applies for all spring rates. Thus, the 
best system performance was obtained with the 
4-inch line, the recirculation valve area at 
twice its nominal value, and the recirculation 
valve spring rate at 1. 2 times its nominal value. 

To determine the consequences of the fail- 
ure of the fuel pump recirculation valve to close 
during startup, the computer model was oper- 
ated with the valve locked open. Since the valve 
is initially open and is closed by fuel pump dis- 
charge pressure overcoming spring forces on 
the valve piston, such a failure might occur as 
a result of a pressure leak at the valve. The 
results show that the engine bootstrapped to the 
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TABLE HI. - EFFECT OF SYSTEM VALVE SPRING 


TABLE II. - RESULTS OF RECIRCULATION RATES ON ACCELERATION TIME 


VALVE FAILURE STUDY 


Engine parameters 

Nominal 

Failure 

p c’ psi 

1040 

788 

P GG’ psi 

1145 

924 

(o/F) TC 

5. 5 

9.4 

(o/F>gg 

0.8 

1.1 


Percent of 

Percent of nominal gas generator 

nominal thrust 

valve spring rate 

chamber valve 
spring rate 

50 

100 

150 


Acceleration time to 90-percent-rated 


thrust chamber pressure, tgg, sec 

50 

— 

1.70 

— 

100 

1.38 

1.60 

2.37 

150 

1.29 

1. 50 

— 


conditions outlined in table n. It is seen that the thrust chamber mixture ratio at these 
conditions is near stoichiometric. Although the computer model had no burnout limits, 
and hence, allowed operation in this region, this occurrence on an engine would result in 
burnout, since this mixture ratio falls outside the design range. 

The sensitivity of the engine start transient to the timing (spring rate) of the remain- 
ing system valves was also determined. The time to 90-percent thrust chamber pressure 
was noted for various gas generator and thrust chamber valve spring rates. The signifi- 
cant results are tabulated in table III. Stall margin was not appreciably affected in this 
study although there was a trend toward longer acceleration times with stiff er gas gener- 
ator valve springs. This was due to the resulting delay in gas generator ignition and 
turbopump spinups. The inverse effect of the thrust chamber valve spring rate is due to 
the increase in flow to the gas generator as the thrust chamber valves are held closed. 

To aid in engine- structure coupling studies (POGO), linear transfer functions to rep- 
resent the M-l system from the pump inlets to the thrust chamber were determined. The 
POGO problem results from the coupling of the vehicle structure and engine dynamics in 
such a way as to produce oscillatory variations in thrust that can be very damaging for 
both manned and unmanned payloads. Determinations of engine transfer functions early 
in the development cycle expedite the solution of the problem since it aids in defining the 
parameters of interest and provides important design data. The transfer functions of in- 
terest were the ratios of thrust chamber pressure to propellant tank pressures, assuming 
sinusoidal pressure perturbations. When the suction lines were removed and the propel- 
lant tanks were essentially placed at the pump inlets, results were obtained that would be 
applicable, regardless of the inlet line dynamics encountered. 

Figures 14 and 15 show the magnitude and phase of the transfer function between 
thrust chamber pressure and oxidizer pump inlet pressure against frequency. Below 
20 cps, the magnitude of the transfer function is approximately constant. Figures 16 and 
and 17 show the magnitude and phase of the liquid oxygen side input impedance against 
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Phase angle, 



Figure 14 - Magnitude of liquid oxygen side transfer func- 
tion between thrust chamber pressure and liquid oxygen 
pump inlet pressure. 



Figure 15. - Phase angle of liquid oxygen side transfer function between 
thrust chamber pressure and liquid oxygen pump inlet pressure. 



Figure 16. - Magnitude of liquid oxygen input impedance 
against frequency for analog-computed and hand- 
calculated cases. 


19 


i 





Figure 17. - Phase angle of liquid oxygen input impedance against frequency 
for analog-computed and hand-calculated cases. 



Figure 1& - Magnitude of fuel side transfer function be- 
tween thrust chamber pressure and fuel pump inlet 
pressure. 



Figure 19. - Phase angle of fuel side transfer function between thrust chamber 
pressure and fuel pump inlet pressure. 
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Figure 20. - Magnitude of fuel input impedance against 
frequency for analog-computed and hand-calculated cases. 



Figure 21. - Phase angle of fuel input impedance against frequency for analog- 
computed and hand-calculated cases. 


frequency. This impedance is defined as the ratio of pump inlet pressure to pump flow 
rate. The effect of cavitation bubbles on the inlet impedance was neglected. As in the 
case of the pressure/pressure transfer function, frequency effects were small below 
20 cps. To support the analog results, hand calculations of the input impedance were 
made and the results of these calculations are also plotted in figures 16 and 17. The 
hand analysis is described in detail in appendix E. 

Figures 18 and 19 show the magnitude and phase of the transfer function between 
thrust chamber pressure and fuel pump inlet pressure. As in the case of the liquid oxy- 
gen side, both computer and hand calculated results were obtained for the fuel flow im- 
pedance Zp and are presented in figures 20 and 21. For frequencies below 50 cps, the 
fuel impedance can be approximated by a constant that is independent of frequency. The 
desired transfer functions between thrust chamber pressure and propellant tank pressures 
can be obtained, for any given suction line impedance, by using the following equations: 




(33) 
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where Z ^ is the impedance of oxidant suction line and 



( 34 ) 


where Zgp is the impedance of fuel suction line. 

Since no combustion dead time was considered in the analysis, it should be noted that 
the addition of an estimated 0. 7 millisecond dead time (ref. 6) to the engine dynamics 
would contribute a maximum of about 30° to the input impedance characteristics at 
100 cps. 


CONCLUSIONS 

The analog computer simulation of the M- 1 rocket engine system yielded the follow- 
ing information concerning engine characteristics: 

1. The effect of ambient pressure on the engine start transient was found to be negli- 
gible with pressures ranging from zero to atmospheric (the figures of merit being fuel 
pump stall margin and the time to reach 90-percent thrust chamber pressure). 

2. Stall margin could, in general, be increased by increasing the recirculation valve 
area at the cost of increased acceleration time. Stall margin was found to be a slight in- 
verse function of recirculation line diameter. A decrease in acceleration time, together 
with an increase in stall margin, was obtained by using the 4-inch line with recirculation 
valve areas and spring rates greater than nominal. 

3. Failure of the recirculation valve to close during the start would result in a steady- 
state thrust chamber mixture ratio near stoichiometric, with a resulting burnout of the 
thrust chamber. 

4. Thrust chamber to propellant tank pressure ratios can be obtained from engine 
impedance characteristics measured on the analog computer. Over the range of ex- 
pected POGO oscillations (0 to 20 cps), the engine impedance characteristics may be as- 
sumed constant with the dominant dynamics expected from the inlet lines. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 11, 1965. 
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APPENDIX A 


SYMBOLS 


A 

cross-sectional area, throat 

K r 

rated turbine static to total dis- 

B 

area, in. ^ 

2 

fluid bulk modulus, lb/in. 

L 

charge pressure ratio 
shaft torque, in. -lb 

C 

2 

fluid capacitance, in. 

J 

hydraulic line inductance, 

9 9 

C D 

orifice discharge coefficient 


sec /in. 

c* 

characteristic exhaust velocity 

l 

length, in. 


(function of chamber mixture 

M 

Mach number 


ratio), in. /sec 

N 

turbopump shaft rotational speed, 

c 

isentropic spouting velocity, 


rpm 


in. /sec 

O/F 

oxidant fuel ratio 

C P 

specific heat of fluid at constant 

P 

2 

absolute total pressure, Ib/in. 

F 

pressure, Btu/(lb)(°R) 
flow correction factor to account 

AP 

absolute total pressure change 

2 

(assumed positive) , lb/in. 


for variations in fluid proper- 
ties 

P 

2 

absolute static pressure, lb/in. 

f i 

general functional relation 

Q 

Q 

volumetric flow rate, in. /sec 


(i = integer) 

q 

rate of heat flow, Btu/sec 

g 

2 

gravitational constant, in. /sec 

R 

universal gas constant, in. /°R 

h 

specific enthalpy, Btu/lb 

R 

hydraulic flow resistance, appro- 

Ah 

specific enthalpy drop (assumed 


priate dimension 


positive) across turbine, 

r 

mean turbine blade radius, in. 


Btu/lb 

S 

Laplace operator, 1/sec 

I 

turbopump moment of inertia, 
2 

T 

total temperature, °R 

J 

in. -lb/sec 

mechanical heat equivalent 
(constant), in. -lb/Btu 

t 

u 

time measured from start, sec 
turbine bucket mean peripheral 

K i 

proportionality constants 
(i = integer), appropriate 

V 

velocity, in. /sec 
plenum volume, in. ® 


dimensions 

W 

rate of fluid weight flow, Ib/sec 
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X general fluid property, appro- 

priate dimension 

y valve stroke, in. 

o 

Z flow impedance, sec/in. 

0 percent burnt products in turbine 

driving gases 

y fluid specific heat ratio 

77 efficiency, (0 < 77 < 1.0) 

p fluid density (constant unless 

otherwise specified), lb/ in. 

r chamber relaxation time, sec 

Subscripts: 

ACT actual valve 

BP gas generator burnt products 

C thrust chamber 

D pump discharge 

d downstream 

F fuel 

FP fuel pump 

FT fuel turbine 

GG gas generator 

HE helium 

ID ideal value 

IL interstage line 

INJ injector 

1 initial value 

J coolant jacket 

J 1 upstream half of coolant jacket 

J 2 downstream half of coolant jacket 

M helium - O/F mixture 


man manifold 

O oxidant 

ORF engine balance orifice 

OT oxidant turbine 

P pump (in general) 

PA piston actuating 

PUV propellant utilization valve 

RC recirculation 

r at rated mixture ratio 

SF fuel suction line 

SK nozzle extension skirt 

SO oxidant suction line 

SV start valve 

T tank 

TAP main line tapoff 

TC thrust chamber 

TUR turbine (general) 

t throat 

u upstream 

V valve 

WG coolant tube wall (combustion 

zone side) 

WL coolant tube wall (coolant flow 

side) 

Wg coolant tube wall (ambient side) 

Superscript: 

- average 
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APPENDIX B 




TIME AND VOLTAGE SCALING 

Since the computers used are 100-volt devices, all system variables must be scaled 
so as not to exceed that value. The voltage representing any variable a is designated 
E a and is related to a by the expression = a/SF where SF is the voltage scale 
factor. To facilitate data recording and the use of servodriven multipliers, the problem 
was slowed down (time- scaled) by a factor of 100. 


Hot Gas and Helium Section 


General parameter, 

a 


Scale factor, Maximum 

SF general parameter, 

a 

max 


^GGO 

0.015 

1 

^GGO 

.8 

80 

W GGF 

1 

100 

W (all others) 

1. 5 

150 

P GG 

12 

1200 

p 

uFT 

12 

1200 

P dFT 

5 

500 

P uOT 

5 

500 

P dOT 

2 

200 

P HE 

50 

5000 

P dSV 

50 

5000 

P d / p u («“> 

.01 

1 

R (all) 

8 

800 

t he 

4 

400 

T (all others) 

40 

4000 

c p (all) 

.04 

4 

( a1i act) f 

8 

800 

( Ah ro) F 

10 

1000 

\ h ACT/ 0 

2 

200 

( Ah ro) 0 

4 

400 
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p 

0.01 

1 

V (all) 

.01 

1 

U F 

30 

3 000 

u O 

5 

500 

C F 

100 

10 000 

c o 

100 

10 000 

N o 

40 

4 000 

n f 

150 

15 000 

l ot 

400 

40 000 

l ft 

400 

40 000 


Main Line and Pump Section 


General parameter, 

a 


Wq (all flows) 
W F (all flows) 

P TO 

P TF 

Pq (all others) 
Pp (all others) 

L OP 

L FP 


Scale factor, 

Maximum 

SF 

general parameter 


a 

max 

40 

4 000 

8 

800 

1 

100 

1 

100 

20 

2 000 

25 

2 500 

400 

40 000 

400 

40 000 


Cooling Jacket and Main Engine 


General parameter, 

Scale factor, 

Maximum 

a 

SF 

general parameter 



a 

max 

T c 

75 

7 500 

t wg 

30 

3 000 

t wl 

30 

3 000 

T W2 

6 

600 

t f 

1 

100 


26 



I 


T J 

2. 5 

2 500 

t inj 

3 

300 

q (all) 

2500 

250 000 

W Q (all) 

40 

4 000 

W F (all) 

8 

800 

w c 

40 

4 000 

P Q (all) 

20 

2 000 

P F (all) 

25 

2 500 

P C 

12. 5 

1 250 
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APPENDIX C 


MIXING OF HELIUM AND GAS GENERATOR BURNT PRODUCTS 


To account for variations in fluid properties due to the mixing of helium and gas gen- 
erator burnt products, it is necessary to derive c , R, and h for the mixture, assuming 
perfect instantaneous mixing in a plenum. The following definitions are used: 


Wj weight of fluid 1 in plenum 

W 2 weight of fluid 2 in plenum 

Wj rate of flow of fluid 1 into plenum 

W 2 rate of flow of fluid 2 into plenum 

Cpi specific heat of fluid 1 

Cp2 specific heat of fluid 2 

Rj gas constant of fluid 1 

R 2 gas constant of fluid 2 

hj specific enthalpy of fluid 1 

hg specific enthalpy of fluid 2 

M subscript denoting desired mixture properties 

X general fluid property (either C , R, h) 

ir 


Then 


*M = 


W 1 X l + W 2 *2 A" 1 ^ + /A dt 

W l +W 2 yCW 1 +W 2 )dt 


Cross multiplying yields 


* u f (W, +W 2 )dt ./>.*!*♦/ 2 dt 

When it is assumed that and X^ do not vary appreciably during the mixing process, 

X uf ,(W 1 +1 * r 2> dt = X lf"l * + *2 f"t dt 
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By differentiating both sides with respect to time, we get 


*i> r i + *2>- VrVi 

Therefore, 

W, W 9 

Xjj, = x, — K- + x _ — K- 

M l W x +W 2 ^W 1 + W 2 

* M ' e*i ♦ (‘ - »* 2 

where /3 is the percent of fluid 1 in the resulting mixture. 

The preceding derivation is based on the assumption that mixing takes place instanta- 
neously in a small plenum. The result may be modified to more closely represent the 
physical process as 


dX 


M 


'M 


dt 


+ * M = /3X 1 +{1 " ^*2 


where is the mixing time constant. 


29 



APPENDIX D 


EQUATION SUMMARY 
Helium and Hot Gas Section 


y GGV 




P dFP dt + *2 


(Dl) 


• / \ r / \ 

“1 1/2 


W GGO = ( C D A ) 0 L 2gp o( AP GGv) 

u 

(D2) 


1 1/2 


W GGF _ ( C D A ) 2gp F( AP GGvJ 


(D3) 

A 0 = f l^ y GGV^ 


(D4) 

A F = V y GGv) 


(D5) 

^ORF " ( C D A ) ORF ^ 2gip ORF AP ORF^ 17 
P GG + P uFT 

(D6) 

/r\n\ 


p ORF 


2R M T M 


(D7) 


_ P uFT A FT ( v P dFT \ ^ r r 
W FT ( K rF p ) 

( r m t m) 7 V uFT/ 


K 0 


'VM 


- V 


+ K„ 






2gy r 


y - 1 

r v 


,(r T .-i)/r t Tl"T 1//2 


1- K. 


rF 


P dFTY' r ~"' r 


P uFT/ 


(D8) 


30 


(D9) 


W, 


P uOT A OT ( k P dOT\ 1//yr 

OT '(%w 1 / 2 V r °w 


K c 




- V 


+ K„ 


y M 




2gy T 


r r - 1 


1- K. 


rO 


P \(y r -i)/y I n > l 1 ^ 2 

*dOT \ r r 


5 uOT/ 


r M 


- 1 


= 8.03x10 


-4 R M 


r M 


'PM 


W 


PuV 


A W 
A PuV w FT 

A FT 


W 


IL 




P dFT “ P uOT " ^IL^Il)^ 


W 


SK 


A SK P dOT / 2 


V(r M -i)/ 2gr X 1 / 2 


( R M T dOT) V2 ^ M + l7 


M 

+ 1 J 


W 


HE 


P HE A SV I 

r p dsv\ 1/, ' HE l 

f 2g y- 

(R T ) 
lK HE A HE> 

We/ 1 

[_ y HE 


AP ORF " P i 


k HE 


j / p dsv\ (yHE ~ 1)/yHE ~ 

. V P HE / 


1 1/2 


( AP GGv) 0 _ P TAPO ■ P GG 


( AP GGv) f " P TAPF " P GG 


GG 


R BP T GG 

V GG 


r 

Jo 


^GGO + W GGF ” W ORF^ dt 


(DIO) 

(Dll) 

(D12) 

(D13) 

(D14) 

(D15) 

(D16) 

(D17) 

(D18) 
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uFT 


^ M f Q ^ORF + ^HE " ^FT " Wp u y) dt 


uFT 


^IL ° 


dFT 


uOT 


,^L £ ^ w OT )dt 


IL 


dOT 


VAm Cv> OT-V 


V dOT ^ 


P dSV _ P uFT + ^SV^HE 


' c pm t m 


1- K 


rF 


P dFT \ 

P uOT/ 


(rM -1 )/^M 


( Ah m) 0 - c pm t uOt 


1- K 


rO 


P 

*dOT \ M 
\iOT/ 


( Ah ACT) F " ^FT^rojp 

( a1i act) 0 = ^0T( Ah ro) 0 

"FT'Mc) I tKs (cL 


"or ' Ks ( c ) 0 + K 6 ( o ) 0 


(D19) 

(D20) 

(D21) 

(D22) 

(D23) 

(D24) 

(D25) 

(D26) 

(D27) 

(D28) 

(D29) 
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7rr OT 

U-O = N 0 
u 30 U 


(D30) 


n rp 

F 30 F 


(D31) 


c O ■ [ 2gJ ( Ah ACT) 0 ] 1//2 

(D32) 

e F = [ 2 gj(Ah ACT ) J 1//2 

(D33) 

L OT = W OT 

( Ah ACT) D 

N o 

(D34) 

T, = _£. W 

( a1i act) f 

(D35) 

FT 2ir FT 

n f 


N ° " ^i~ fo (L ° t ‘ L op )dt 


(D36) 


N, 


= ^Lfo ( L FT- L Fp) dt 


(D37) 


T uOT _ T M 


"PM 


^FT^ACt), 


T dOT ~ T uOT 


W FT + W PuV 


( a1i act). 


"PM 


(D38) 


(D39) 


0 = W GGO + W GGF 


W GGO + W GGF + W HE 


(D40) 
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C PM ~ /3c pBP + ^ " ^ c pHE 

(D41) 

% = ^^BP + ” ^ R HE 

(D42) 

c PM T M = i3c pBP T BP + ^ " ^ c pHE T HE 

(D43) 

T p 1-y HE //y HE T p 1_r HE/ r HE 
1 HE J ^HE “ ii 

(D44) 


Main Line and Pump Section 


w 


so 


— r 

J so Jo 


p _ p - idt 

p TO F SO „ w SO) ai 


(D45) 


W 


SF 


— f 

•^SF Jo 


1 ( P TF - P SF - "If )* 


(D46) 


W 


OP 


y ~ y*i 

z SO -/o 


P on " P + A P/ 


SO IT t2 


Idt 


(D47) 


W 


FP 


— c 

^SF - A ) 


P -P + AP - ^dt 

P SF P dFP +A± FP o W FP/ ar 


(D48) 


W 


DO 




P dOP “ P TAPO " ^DO^ DO/ dt 


(D49) 


W 


DF 


.- i _ r 

•^DF *4) 


P dFP “ P TAPF ^DF^Df)^ 


(D50) 


/ . \ / \ 

/ \ 

( W TCv) 0 = ( C D A TCV/ 0 

2 gP 0 ( AP TCv) 0 


1/2 


(D51) 
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(w TCV ) - ( C D A TCvL 2gp F( AP TCv) 


Wrcv = (c d a ) rcv ( 2 ^f ap rcv) 1/2 


A TCVO “ f 3^ P dOP + P dFp) 


V TCVF “ f 4^ P dOP + P dFP^ 


^RCV “ f 5^ P dFP^ 




unstalled 




= K 14 + K 15 


stalled 


K 16 + K 17 


stalled 




iP op _ y ,r r opV, w op 
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/ ^o-wopMt 

(D63) 

P SF ' Pf v sf f 0 ^F W FP> dt 

(D64) 

PdOP %o ^ ^ " WDO>d ‘ 

(D65) 

PdFp -p F v DF X (Wpp WDF>dt 

(D66) 

P TAPO = PqVdo f 0 [ W DO“ W GGO ( W TCv) 0 _ 

dt (D67) 

P TAPF v L [ W DF " W GGF " ( W TCv) f 

Pf v df 0 L * J 

dt (D68) 

*uJ = P J + -f W Jl 

(D69) 

( AP TCv) g = P TAPO “ ( P INj) 0 

(D70) 

( AP TCv) f = P TAPF " P uJ 

(D71) 

AP RCV = P dFP P TF 

(D72) 

Cooling Jacket and Main Engine (see fig. 7) 


9l - K 24< W C>°’ 8 < T C - T WG> 

(D73) 

q 2 = ^S^WG “ T WL^ 

(D74) 


q 3 “ ^e^TCv)^. 8 ( T WL “ T J^ 

(D75) 

q 4 = K^y^tcv)^' 8 ^ T W2 ■ T j) 

(D76) 

T WG = *^28 f 0 ( q l ' q 2^ dt 

(D77) 

T WL = ^9 ^ q 2 " q 3^ dt 

(D78) 

T W2 = K 30 f Q q 4 dt 

(D79) 

— q-j + q 4 

T T = Ko i / . ’ +T-, 

J 31 ( w tcv) f F 

(D80) 

( t inj) f - 2 *J - t f 

(D81) 

^ W INj) 0 = ( W TCv) 0 

(D82) 

W J1 = ( w tcv) f 

(D83) 

^J2 R j**J ( P INj) F 

(D84) 
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APPENDIX E 


CALCULATION OF M-l OXIDIZER FLOW IMPEDANCE 


To aid in engine-vehicle coupling studies (POGO), it is necessary to obtain thrust to 
propellant tank pressure frequency response data. A primary step in obtaining this in- 
formation is the calculation of flow impedances (Z Q = 9P u op/ 9W op) * or the various seg- 
ments of the M-l system. The segment considered in this analysis is illustrated in fig- 
ure 22. When standard perturbation techniques and the design values for flow rates and 
pressures are used, the electrical analog to the fluid system shown in figure 23 can be 
obtained. One tt- shaped section was used to approximate the discharge line. All pres- 
sures and flows represent perturbations from the operating point. The values computed 
for the system parameters are shown in table IV. These values, together with classical 
circuit reduction techniques, allowed us to obtain the following expression for the oxidant 
impedance Zq: 


where 


z o = 


K 33 s4 + K 34 s3 + k 35 s2 + K 36 S + *37 
K 00 S 4 + K on S 3 + K.nS 2 + K„ i S + K 


38 L 


39 k 


4<T 


4P 


42 


K 33 = 7. lx 10" 13 

K 37 = 2. 92 

K 4q = 1. 17X10" 5 

K 34 = 6.79X10" 9 

K 38 = 2. 85X10“ 12 

K 41 = 7.35X10" 3 

K 35 = 1.64X10' 5 

Kgg = 1. 51X10" 8 

K 42 = 3. 43 

K 36 = 1. 42X10" 2 





( W|NJ>F 


Figure 22. - Block diagram of M-l oxidizer system. 
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TABLE IV. - SYSTEM PARAMETER VALUES 


System parameters 


Computed 

values 



Definition 

Equation 


2 

Liquid oxygen pump flow resistance, sec/in. 

* _ s A? OP 

op aw 

ow 0 p 

0.25 

Liquid oxygen discharge line flow resistance, 
sec/in. 2 

g _ sap do 

DX 0W 

dW DO 

1.15X10" 2 

Liquid oxygen thrust chamber valve flow 
resistance, sec/in. 2 

„ S ap tcvo 

xCVO 

dw TCVO 

3. 54xl0" 2 

2 

Thrust chamber flow resistance, sec/in. 

«C=". c 

3W TCVO 

0.292 

Liquid oxygen injector flow resistance, 
sec/in. 2 

D 9 ap injo 

IN JO 9^ 

0W TCVO 

0.276 

Liquid oxygen discharge line inductance, 

2/. 2 
sec /in. 

, 1 DO 

DO " A e 
A DCr 

2.96X10' 3 

2 

Fluid capacitance of discharge line, in. 

0 P 0 V D0 

DO R 

B o 

1. 9xlO' 3 

2 

Fluid capacitance of thrust chamber, in. 

c = — 
R c 

3. 43X10' 3 


A similar procedure was followed to determine the impedance for the fuel system when 
heat transfer in the jacket was neglected. The resultant flow impedance was 

z _ K 43^ + K 44 s5 4 K 4 5 s4 + + k 47^ + K 48 S * K 49 

F K 50 s6 + K 51 s5 + K 52 s4 + K 53 s3 + K 54 s2 + K 55 S + K 56 

where 


K 43 = 3. 35x10" 17 

K 48 = 4. 66X10“ 1 

K 53 = 1. 99X10 

K 44 = 3. 78X10" 13 

K 49 = 25.4 

K 54 = 9.72x10 

K 45 = 1. 06x10“ 9 

K 50 = 5. 7 IX 10" 18 

K 55 = 7. 39X10 

K 46 = 1.38X10" 6 

K 51 = 6. 32x10" 14 

K 56 = 3. 52 

K 47 = 7. 94X10" 4 

K 52 = 1.66x10" 10 
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APPENDIX F 


M-l START AND SHUTDOWN SEQUENCES 
Start Sequence 

At arm switch on, the propellant utilization valve is at 22. 2 percent of its steady- 
state open area. The fuel pump circulation valve (recirculation valve) is open (spring 
loaded). Automatic sequencing controls control the engine startup as follows: 

(1) At fire switch one, the start valve is opened, and the gas generator valve pilot 
valve, thrust chamber, and gas generator igniters are energized. The propellant utiliza- 
tion valve starts moving toward its steady- state position. 

(2) High pressure helium gas is admitted to the turbines initiating pump rotation. 

(3) The thrust chamber valve opening is initiated by pump discharge pressures. 

(4) The gas generator valves are actuated by increasing hydrogen pressure taken 
from the gas generator propellant line. A timing orifice controls the opening rate. 

(5) Propellants enter the gas generator combustion chamber and are ignited. 

(6) Propellants enter the thrust chamber combustion chamber and are ignited. 

(7) Rising fuel pump discharge pressure closes the spring loaded fuel pump circula- 
tion valve. 

(8) As gas generator chamber pressure increases, the helium check valve is closed, 
which terminates the helium flow from the start bottle to the turbines. 

(9) The start valve is closed at 90-percent-rated thrust chamber pressure, which 
vents the helium start lines. 

(10) The gas generator and thrust chamber igniters are deenergized at 90- percent- 
rated thrust chamber pressure. 

(11) The engine bootstraps to its rated operating conditions. 


Shutdown Sequence 

After an initiating input signal, automatic sequencing controls control the shutdown 
as follows: 

(1) The shutdown signal deenergizes the gas generator valve pilot valve to the vent 
position, which vents the actuator and allows spring forces to close the valves. 

(2) The closing of the gas generator valves stops the flow of propellants to the gas 
generator and, consequently, the supplying of drive gases to the pump turbines. 

(3) With a loss of driving power, the pumps begin decelerating and the discharge 
pressures decay. 
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(4) Decaying pump pressures allow spring forces to close the thrust chamber valves 
and open the recirculation valve. 

(5) As the thrust chamber and gas generator valves close and propellant flow is shut- 
off, the turbopumps stop and engine operation is terminated. 
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